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We study a deterministic scale-free network recently proposed by Barabasi, Ravasz and Vicsek. 
We find that there are two types of nodes: the hub and rim nodes, which form a bipartite structure 
of the network. We first derive the exact numbers P(k) of nodes with degree k for the hub and rim 
nodes in each generation of the network, respectively. Using this, we obtain the exact exponents 
of the distribution function P(k) of nodes with k degree in the asymptotic limit of k — > oo. We 
show that the degree distribution for the hub nodes exhibits the scale-free nature, P(k) oc A; -7 with 
7 = In 3/ In 2 = 1.584962, while the degree distribution for the rim nodes is given by P(k) oc e ' k 
with 7' = ln(3/2) = 0.405465. Second, we analytically calculate the second order average degree 
of nodes, d. Third, we numerically as well as analytically calculate the spectra of the adjacency 
matrix A for representing topology of the network. We also analytically obtain the exact number of 
degeneracy at each eigenvalue in the network. The density of states (i.e., the distribution function 
of eigenvalues) exhibits the fractal nature with respect to the degeneracy. Fourth, we study the 
mathematical structure of the determinant of the eigenequation for the adjacency matrix. Fifth, we 
study hidden symmetry, zero modes and its index theorem in the deterministic scale-free network. 
Finally, we study the nature of the maximum eigenvalue in the spectrum of the deterministic scale- 
free network. We will prove several theorems for it, using some mathematical theorems. Thus, we 
show that most of all important quantities in the network theory can be analytically obtained in 
the deterministic scale-free network model of Barabasi, Ravasz and Vicsek. Therefore, we may call 
this network model the exactly solvable scale-free network. 

PACS numbers: 89.75.-k,89.75.Da,05.10.-a 



I. INTRODUCTION 

There has been a notable progress in the study of the 
so-called scale-free network (SFN) 0,011110,11111111 
Hcj for recent years. In the network theory, the random 
network model was first invented by Erdos and Renvipd|. 
Recently it was generalized to the small world network 
models [H, El, Q EE EE E3, El E3, Further- 
more, about five years ago, the SFN was discovered by 
studying the network geometry of the internet [lL 12.1a. l2ll 
I2I I23I. I24L Faloutsos brothers and Albert, Jeong and 
Barabasijl |H H |H ^ first showed the scale-free 
nature of the internet geometry and opened up an area 
for studying very complex and growing network systems 
such as internet, biological evolution, metabolic reaction, 
epidemic disease, human sexual relationship, and econ- 
omy. These are nicely summarized in the reviews by 
Barabasi0. 

As was studied in the literature^, the nature of these 
SFNs is characterized by the power-law behavior of the 
distribution function. Here the number of nodes with 
order k can be fit by 



P(k) oc fc~ 7 , 



(1) 



where 7 ~ 1 — 4. In order to show the power-law dis- 
tribution of the SFN, Albert and Barabasi first proposed 
a very simple model called the Albert-Barabasi (AB)'s 
SFN model0 EJ |H |H || . 

This system is constructed by the following process: 
Initially we put mo nodes as seeds for the system. Every 



time a new node is added m new links are distributed 
from the node to the existed nodes in the system with a 
preferential attachment probability 



(2) 



where fcj is the number of links at the i-th node. The 
development of this model is described by a continuum 
model 



dki tt / j n mki 

= mil, ki = . 

dr V ' 2t 



(3) 



Then at time r the system consists of N(r) nodes and 
the L(t) links with 



L{t) 



i5> 



(4) 



i=l 



As studied by Barabasi et. al.0 this model exhibits an 
exact exponent of 7 = 3 for the power-law. Thus, it 
has been concluded that the essential points of why a 
network grows to a SFN are attributed to the growth of 
the system and the preferential attachment of new nodes 
to old nodes existed already in the network. 

From the above context, time evolution to construct 
an SFN has been intensively studied in the ^4i?-model as 
well as other models. And many works have appeared, 
regarding nodes and links as metaphysical objects such 
as agents and relationships in an area of science |2|. How- 
ever, most approaches were based on the numerical ap- 
proach. And the spectra of the adjacency matrix A for 



2 



the SFN have been studied numerically S IM |H E3 • 
Therefore, apart from the purposes for the numerical 
analysis, the continuous-time SFN models such as the 
AS-modcl arc not good enough to see what is going on 
in the network geometry in the microscopic level. 

Instead of such a continuous-time SFN model, a new 
type of the SFN models, sometimes called determin- 
istic or hierarchical SFN models, has been proposed 
by Barabasi, Ravasz, and Vicsek[28| and Ravasz and 
Barabasi^ (We would like to call it the DSFN in this 
paper). In the former, the study showed a power-law be- 
havior of the network analytically, while in the latter, the 
study showed it numerically. However, much is still lack 
and left unanswered. 

On the other hand, there is a very important problem 
on the maximum eigenvalue, X ma x, of the adjacency ma- 
trix, A. As was numerically studied |25l 129. 127| . the max- 
imum eigenvalue for the AB-model is bounded by \Jk max 
such that A 

max — const .^J 7c maa; , where fc ma;E means the 
maximum order of nodes. And the numerical studies 
showed that k ma x °c \/~N. And therefore, the numerical 
studies validated \ ma x < const. N 1 / 4 . On this problem, 
very recently, Chung, Lu and Vu|3l| have proved a very 
general theorem: In a complex network model, X m ax is 
always bounded by the lower and upper bounds such that 
a < Xmax < b. Define the second order average degree of 
nodes, d (see Sec. IV). Then, (CI) if the exponent 7 > 2.5, 
then const.d < X max < const.y/k max . (C2) If the expo- 
nent 7 < 2.5, then const. \Jk max < X max < const.d. (C3) 
And if the exponent 7 = 2.5, then a transition happens. 
From this theorem, the AB-model belongs to the first 
category since 7 = 3. However, in spite of the seemingly 
important theorem, because of the lack of other good ex- 
amples other than the AB-model, the examples for the 
other categories have not yet been so well-known so far. 

So, the purpose of this paper is to study in much de- 
tails the DSFN model proposed by Barabasi, Ravasz, and 
Vicsek 28] in order to give a good example for the other 
category of the theorem. This study will provide a rigor- 
ous treatment for the complex network model. 

The organization of the paper is the following. In Sec. 
II, we will introduce the DSFN model that was first stud- 
ied by Barabasi, Ravasz, and Vicsek 28]. In Sec. Ill, 
we will present the exact numbers of nodes and its de- 
grees. And we calculate the exact scaling behavior of the 
nodes. In Sec. IV, we will calculate the exact number of 
the second order average degree d, using the exact num- 
ber distribution of nodes and degrees. In Sec. V, we 
will introduce the so-called adjacency matrix in the net- 
work theory to the DSFN. In Sec. VI, we will introduce 
the eigenvalue problem of the adjacency matrix. In Sec. 
VII, we will present the numerical results of the spectra 
of the adjacency matrices for the DSFNs up to the n = 7 
generation. And we will derive the exact numbers of de- 
generacies in the spectrum. In sec. VIII, we will present 
an analytical method that deduces the exact numbers of 
the degeneracies in terms of the irreducible polynomials 
for the DSFN. We will also present some conjectures for 



the polynomials. In Sec. IX, we will discuss the role of 
the roots of the irreducible polynomials. In Sec. X, we 
show that there is a hidden symmetry in the adjacency 
matrix of the DSFN. In Sec. XI, we will study the zero 
modes and its index theorem in the DSFN. In Sec. XII, 
we will discuss the nature of the maximum eigenvalue in 
the spectrum of the DSFN. We will prove several theo- 
rems using some mathematical theorems. And we will 
discuss the relationship between the results of our theory 
and the Chung, Lu, and Vu's theorem|3l|. In Sec. XIII, 
conclusion will be made. 



II. DETERMINISTIC SCALE-FREE NETWORK 

Let us introduce the DSFN model invented by 
Barabasi, Ravasz, and Vicsek 28]. The development of 
this network is shown in Fig.l. The black and red nodes 
show the hub and rim nodes. We call the most connected 
hub and rim the root and leaf, respectively. 
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FIG. 1: (Color online) The deterministic scale-free network. 
The black and red nodes show the hub and rim nodes. We 
call the most connected hub and rim the root and leaf, re- 
spectively. This network is a bipartite structure. 

In this model the total number of nodes, N(n), the 
total number of links, L(n), and the maximum number 



3 



of links, k max (n) are given by 

N(n) = 3", (5) 

L(n) = 3L(n - 1) + 2" = 2(3" - 2"), (6) 

k max (n) = 2 + 2 2 + ■ ■ • + 2" - 2(2" - 1). (7) 

respectively. From these we find the their development 
as follows: As n = 0, 1, 2, 3, 4, • • • , then 



AT(ra) = l,3,9,27,81 r 



L(n) = 0,2,10,38,130,' 



U(n) =0,2,6,14,30,- 



(8) 

(9) 
(10) 



Let us consider the average link number (i.e., the av- 
erage degree) (k) of a network. It is defined by 



N(r 



(*,)(„) = -L-^k 

N(n) 



(11) 



The meaning of this is just the number of links per node 
(i.e., the average degree). We may call it the first order 
average degree. On the other hand, if we use the num- 
ber P(k) of nodes with k degree, then we can write the 
average as 



k maa: (n) 



(12) 



fe=i 



Therefore, we find that the conversion is carried out by 

(13) 



N{n) 

P(k) = £ s kM 

i=l 



Now, we are able to calculate (k) for our DSFN. Substi- 
tuting Eq.s (8) and (9) into Eq.(ll), we obtain 

N(n) 



2L(n) 

N(n) ^ "* ~ mn) 



2 • 2(3" - 2") 
3" 



= 4 1 



4. 



(14) 



In this way, even though the network becomes very 
complex as n — > oo, the average approaches a finite con- 
stant 4. This is due to the following fact: In this DSFN 
as the iteration is repeated, the order of the most con- 
nected hub becomes large indefinitely while its number 
remains very few (i.e., « o(l)). On the other hand, the 
numbers of the very few connected nodes become large 
indefinitely. Hence the sum of the magnitude of order (k) 
times the number (P(k)) of nodes with the order k may 
remain finite. 



III. EXACT NUMBERS OF NODES AND 
DEGREES 

Let us find the exact numbers of nodes and degrees. 
This would be very crucial for our later purposes in order 
to evaluate many quantities in the network theory. 

As was discussed by Barabasi, Ravasz, and Vicsekpsj. 
in the DSFN there are two categories of nodes called 
"hub" nodes and the "rim" nodes. As they called the 
most connected hub node the "root" (shown as black 
dots in Fig.l), we would like to call the most connected 
rim node the "leaf (shown as red dots in Fig.l). From 
seeing Fig.l, the locations of the root node and the leaf 
nodes look very similar to those of a hub and rims in an 
umbrella. While there exists only one roof node in each 
generation of the network, the number of the leaves can 
increase very rapidly. 

Let us first consider the hub nodes. In the i-th step, 
the degree of the root, is 2' 1+1 — 2. In the next iteration 
two copies of this hub will appear in the two newly added 
units. As we iterate further, in the n-th step 3" _l copies 
of this hub will be present in the network. However, the 
two newly created copies will not increase their degree af- 
ter further iterations. Therefore, after n iterations there 
are 2 • "y 1 ^ 1 ^ 1 nodes with degree 2 i+1 — 2. 

Let us next consider the rim nodes. In the i-th step, 
the degree of the most connected rim, the leaf, is i. And 
the number of the such nodes is 2\ In the next iteration 
one copy of the leaves will be kept the same and two 
copies of the leaves will appear in the two newly added. 
As I iterate further, in the n-th step 3" _I copies of the 
leaves will be present in the network. Therefore, after n 
iterations there are 2 l • 3™- 2 - 1 nodes with degree i. 

Now, denote by k the degree of nodes and denote by 
P(k) the total number of nodes with degree k. Hence 
we get the following: 



root nodes 


leaf nodes 


k 


P{k) 


k 


P(k) 


2 


2-3 71-2 


1 


2 • 3"^ 


6 


2 • 3"~ 3 


2 


2 2 • 3"~ 3 


14 


2 • 3"~ 4 


3 


2 A ■ 3"~ 4 










2 l+1 - 2 


2 • 3™" 1 -* 


i 












2 ,w - 2 


2-3 


n-2 


2 n -' 2 ■ 3 


2" - 2 


2 


n- 1 


2 n-l . ! 


2 «+i _ 2 


1 


n 


2" 



TABLE 1. The number P(k) of nodes with degree 
k for the root nodes and leaf nodes. 

As was shown by Barabasi, Ravasz, and Vicsek[28j. 
consideration for the root nodes is enough to derive the 
scaling exponent of the distribution function P(k) for the 
root nodes in the network. Picking up the 2-3" -1- * nodes 



4 



with degree 2 l+1 — 2, we can regard P (fc) as 2-3™ 1 1 
and k as 2 i+1 — 2. Then eliminating i, we can derive 



where 



P{k) oc fc~ 7 , 



T = --| = 1.584962. 
In 2 



(15) 



(16) 



This shows a scale-free nature (i.e., the fractal nature) of 
the hub nodes in the network as we expect |30|. 

On the other hand, it is not true for the leaves in the 
network. In this case, regarding P(k) as 2 l ■ 3"- 1 -* and 
k as i, we find 



/>(/,.) x [| ) = e^' k 



where 



7 = In ( - ) = 0.405465. 



(17) 



(18) 



This shows that the scaling nature of the leaf nodes is 
not scale-free but exponential. 

In this way, the scaling nature of the roots and that 
of the leaves in the DSFN are different from each other. 
Thus, we are led to a a certain model which consists of a 
multifractal nature of the complex networks. 



IV. THE SECOND ORDER AVERAGE DEGREE 

Let us calculate the second order average degree, d(n). 
It is defined by 

N{n) AWxO) 

L(n) ^ L(n) ^ 

Th is q uantity was recently introduced by Chung, Lu and 
Vu|3l|]. Roughly speaking, the meaning of this quality 
is the average degree per link. In other words, it is the 
average degree weighted with the preferential attachment 
such that 



N(n) 

d(n) = 2J kiUi(ki). 



(20) 



From Eq.(14) we have 



N(n) 

<*»<»>- 555 E' 



2L(n) 



we derive 



d(n) 



, (fc 2 )(n) 
' (k)(n) ' 



(21) 



(22) 



where (k 2 )(n) is defined by 



(k 2 )(n) 



N(n) 



N(n) 



,(„) 



i=l 



N(n) 



E fc2p ( fc )' ( 23 ) 



fc=i 



the second moment per node. As was shown before, the 
average degree (k)(n) converges to 4 as n — > oo. Hence, 
the second order average degree d(n) becomes propor- 
tional to the second moment (k 2 )(n) in the limit such 
that 



(fc 2 )(n) 



(24) 



Before going to calculate the second order average de- 
gree, let us first check whether or not the distributions 
given in the tables reproduce the correct results for the 
total numbers of nodes and links. This problem is trivial. 
But as we will see, this is very instructive for our purpose 
here. 

Let us show that the distributions of the nodes and 
degrees are exact. Let us first calculate the total number 
of nodes, once again. To do so, we just sum them up as 
follows. For the root nodes we obtain the following: 



Y P(k) = l + 2 + 2-3 + --- + 2-3 



n-2 



k^root 



= 1 + 2(1 + 3 



) = 1 + 2 



in— 1 



- 1 



3-1 



= 3 



n-l 



(25) 

This is the total number of the root (hub) nodes (i.e., the 
total number of black dots in Fig.l). For the leaf nodes 
we have the following: 



P{k) = 2-3"" 2 +2 2 -3 

k(^leaf 

= 3"- 1 +2-3"- 2 +2 2 -3 n ~ 3 



n— 3 



in— 2 



3+2"- 1 - 1+2" 



-2"- 2 -3+2 n - 1 -l+2"-3"- 1 



3" - 2" 



- 2" - 3"- 1 = 3" - 3"" 1 = 2 • 3"- 1 . 



(26) 



This is the total number of the leaf (rim) nodes (i.e., 
the total number of red dots in Fig.l). Hence, by the 
addition of both we obtain the desired result: 



J2 p (k) = 3" = N(r 



(27) 



This proves Eq.(5). 

Let us next consider the total number of links in the 
network. We calculate it for the root and leaf nodes 
separately. For the root nodes, we have the following: 

J2 kP(k) = (2" +1 - 2) • 1 + (2 n - 2) • 2 + (2™- 1 - 2) • 2 • 3 

k£root 



■5 



+ • • • + (2 s - 2) • 2 • 3™~ 3 + (2 2 - 2) • 2 • 3 



371 — 2 



+x"~ 2 • 3 + x"- 1 • l + x" 



2"+ 1 .i + 2".2 + 2"- 1 -2-3 + - • - + 2 3 -2-3"- 3 + 2 2 -2-3 r ' 



a;(1.3«-^ +a; .3«--i + ^.3n-4 + . . . +x n-J. 3+x n-2^ +x 



- (2 + 2 2 + 2 2 • 3 H h 2 2 • 3"~ 3 + 2 2 • 3™~ 2 ) 



2«+i _ 2 _|_ 2 3 (2™ 2 + 3 • 2" ' 5 + • • • + 3" 3 • 2 + 3" 2 ) By differentiating this with respect to x, we obtain 



->n— 3 



jn—3 



in-2 



3 - X 



(33) 



-2 2 (l + 3 + --- + 3 



n-2\ 



/'(*) 



3- x 



2 n+l _ 2 + 2 3 



3 n-l _ 2 n-l 

3-2 



gn-l _ X 

3-1 



+X 



-(n- l)x r > 



-2 2t n ^ x n ^ 



3 — x 



(3-x) 2 



= 2(3" - 2"). (28) 
Similarly, for the leaf nodes, we have the following: 

kP(k) = 1 • 2 • 3™~ 2 + 2 • 2 2 • 3"~ 3 + 3 • 2 3 • 3"~ 4 

k£leaf 



+nx n -\ (34) 
Therefore, substituting x = 2 into Eq.(34), we obtain 

3 n— 1 2 n — ^ 



/'(2) 



3-2 



+ --- + (n- 2)2™~ 2 ■ 3 + (n - 1)2™- 1 • 1 + n2 n 



= 2(1 • 3"~ 2 + 2 • 2 • 3"~ 3 + 3 • 2 2 • 3™~ 4 + • 



+(n - 2)2™- 3 ■ 3 + (n - 1)2™- 2 • 1 + n2 n - 1 ). (29) 

At first glance, the summation of Eq.(29) seems quite 
tough. But to perform this calculation we can use a 
mathematical trick as follows: Define a function f(x) 
of x such that 

n — 2 | 2 on— 3 i _3 on — 4 



f{x) = x ■ T- 2 + x 2 • 3™~ 3 + x 3 • 3™~ 4 + 



+x n - 2 -3 + x™- 1 -l + x™. 



Differentiating this with respect to x, we obtain 
f'(x) = 3"~ 2 + 2x • 3"~ 3 + 3x 2 • 3"~ 4 + • • 



(30) 



+(n - 2)x™- 3 ■ 3 + (n - l)x n - 2 -1 + nx 



„ra-l 



(31) 



Substituting x = 2 into Eq.(31) and comparing it with 
Eq.(29), we find 



£ fcP(fc) = 2/ (2). 



(32) 



Therefore, once /(x) is found as a simple function of x, 
we can calculate the sum. This can be done as follows. 

f(x) = x • 3™~ 2 + x 2 • 3™~ 3 + x 3 • 3™~ 4 + • • • 



+2 



-{n- 1)2™- 2 3"" 1 - 2™- 1 



3-2 



(3 - 2) 2 



+n2 



n-l 



= (3™- 1 -2"- 1 )+2 [-(n - 1)2™- 2 + 3"- 1 - 2™- 1 ] +n2 n - x 

= 3" - 2". (35) 
Substituting this into Eq.(32), we obtain 



J2 kP(k) = 2(3" -2"). 



(36) 



k£leaf 



In this way, explicitly using the exact numbers of nodes 
and degrees, we can show that each sum produces the to- 
tal number of links in the network as we expected. Hence, 
this proves Eq.(6). This situation encourages us to per- 
form the calculation of the second order average degree 
of Eq.(19). 

Let us next do this. In Eq.(19), we need to separate it 
into two parts of the sum as follows: 

JV(n) k max (n) 
x ; i— 1 



Jfc=l 



L(n) 



k 2 p(k)+ k2p ( k ) 



(37) 



K k£root 



kGleaf 



Let us consider the sum for the root nodes. Using the 
previous result, we obtain the following: 

k 2 P(k) = (2 n+1 -2) 2 -l+(2"-2) 2 -2+(2"- 1 -2) 2 -2-3 

kEroot 



H h (2 3 - 2) 2 • 2 • 3"~ 3 + (2 2 - 2) 2 • 2 • 3™~ 2 



(2™ +1 - 2) 2 



+2[(2 2 )™ + (2 2 )™- 1 3 + • • • + (2 2 ) 2 3™- 2 ] 



8 [2« + 2"" 1 3 + • • • + 2 2 3"" 2 ] 



+8[l + 3 + --- + 3™~ 2 ] 



(2™ +1 - 2) 2 



+2 5 [(2 2 ) n " 2 + (2 2 )™~ 3 3 + • • • + 2 2 3™~ 3 + 3"~ 2 ] 



+(n - 2)x™- 2 ■ 3 + (n - l)^" 1 • 1 + nx n . 
Differentiating this with respect to x yields 

(xf'(x)Y = 3"~ 2 + 2 2 x • 3 n ~ 3 + 3 2 x 2 • 3 n ~ 4 + • • 

+(n - 2) 2 x"- 3 ■ 3 + (n - l) 2 x"- 2 • 1 + nV -1 . 
Therefore, substituting x = 2 into this, we get 

| x=2 = 3™~ 2 + 2 2 • 2 • 3"" 3 + 3 2 • 2 2 • 3™~ 4 4 

+(n - 2) 2 2"~ 3 ■ 3 + (n - 1) 2 2"- 2 • f + nV -1 . 
Comparing this with Eq.(39), we get the relation 

k 2 P(k) = 2(xf(x)Y | x=2 . 



Let us evaluate this. Using Eq.(34), we obtain 

xf'(x) = X 



3 n-l _ jn-l 



3 -a; 



2 5 [2™~ 2 + 2"~ 3 • 3 H h 2 • 3"~ 3 + 3™~ 21 



+2 3 [l + 3 + --- + 3™" 2 ] 



(2" +1 - 2) 2 + 2 5 



(2 



2\n — 1 on — 1 



2 2 -3 



3«-i_ 2 "- 1 ,3 n_1 -l 
-2 5 h 2 3 



3-2 



3-1 



: (2 rl+1 - 2) 2 + 2 5 (4™~ 1 - 2 • 3™ _1 + 2" _1 ) + 2 2 (3"~ 1 - 1) 

(38) 

Let us consider the sum for the leaf nodes. 

k^leaf 



+x' 



-{n-l)x r ' 
3 — x 



+ 



(3-z) 2 



Differentiating this with respect to x, we have 



in- 1 ~n — 1 



on— a. _ 

(*/'(*))' = ' 



+3x 



-(n — l)x 



3-x 



n — 2 on — 1 sr,n— 1 



3 — x 



(3-x) 2 



-(n - l)(n - 2)x™- 3 -2(n - l)x"- 2 



3 — x 



(3-x) 2 



= l 2 • 2 • 3™ -2 + 2 2 • 2 2 • 3™~ 3 + 3 2 • 2 3 • 3"~ 4 



+(n - 2) 2 2"- 2 ■ 3 + (n - 1 ) 2 2™- 1 • f + n 2 2" 



= 2(1 2 • 3™~ 2 + 2 2 • 2 • 3™~ 3 + 3 2 • 2 2 • 3™~ 4 + 



+x' 



+ n 2 x n ~ 1 . 



(3 - x) 3 

Substituting x = 2 into this, we get 

^n— 1 2 n— 1 



(*/'(*))' |x=2 = 



3-2 



+(n - 2) 2 2"~ 3 ■ 3 + (n - 1 ) 2 2™- 2 • f + n z 2 n ~'). (39) 

To calculate Eq.(39) we need a similar trick as before. 
Using Eq. (31), we can define 

xf(x) = xT- 2 + 2x 2 • 3"" 3 + 3x 3 • 3"~ 4 + ■ ■ ■ 



+3-2 



-(n- 1)2™- 2 3"" 1 - 2™- 1 



3-2 



(3-2) 2 



+2 2 



-{n- l)(n-2)2"- 3 -2(n - 1)2 



3-2 



+ 



(3-2)2 



7 



)TO— 1 " 



(3-2)3 

= 3»-i _ 2™" 1 + 6[-(n - 1)2"" 2 + - 
-2 2 [-(n- l)(n - 2)2™~ 3 - 2(n - 1)2™~ 2 + 2(3 n_1 - 2"" 1 )] 



V2"" 1 



_ 3«-i _ 2™" 1 - 6(n - l)2 n_1 + 6(3 n_1 - 2"^ 1 ) 

-(n-l)(n-2)2 Tl - 1 -2(n-l)2"+2 3 (3 Tl - 1 -2"- 1 )+n 2 2™- 1 

= (1 + 2 + 2 2 + 2 3 )(3™" 1 - 2™- 1 ) - 3(n - 1)2™ 

+{n 2 - (n - l)(n - 2) - 4(n - 1)}2™~ 1 

= 5 • 3" - (2n + 5)2™. (46) 
Substituting this into Eq.(43), we get 

5^ fc 2 P(fc) = 2(5 ■ 3™ - (2n + 5)2"). (47) 

Using Eq.(47) together with Eq.(38), we obtain 

k 2 P(k) = (2™ +1 - 2) 2 + 2 5 (4™~ 1 - 2 • 3™- 1 + 2™" 1 ) 

k£all 



+2 2 (3 n_1 - 1) + 2(5 • 3" - (2n + 5)2"). (48) 

Therefore, the above procedure enables us to evaluate 
the second order average degree as follows: 



1 



->n+l i-y\2 | ( ATI— 1 _ 2 . 3 n— ^ _j_ 2 n_ ^ 



4(3" - 2") P 2) + 2 (4 " 

+2 2 (3"- 1 - 1) + 2(5 ■ 3™ - (2n + 5)2™)] 



1 



2/ V3 



"5 



(49) 



In this way, the second order average degree d is calcu- 
lated explicitly. We are able to show that it diverges as 
an exponetial-law. 



V. ADJACENCY MATRIX 

Let us consider the adjacency matrix, A, in the theory 
of network [ll I2I I2H I29I |33J | . This matrix is very important 
for the theory of network, since it represents the topology 
of the network structure. Denote by a,j the component 
corresponding to the link between the z-th node and the 
j-th node. In the theory of network, the elements have 
only or 1 according to whether or not there is a link. 
Therefore, in general, the adjacency matrix is defined by 



1 if there is a link 
otherwise 



(50) 



Let us consider the adjacency matrix for the DSFN. 
Denote by A n the adjacency matrix for the n-th genera- 
tion of the network. Using the numbering given in Fig.l, 
the adjacency matrix is defined by the following. 



(51) 



Ao 












1 


1 












( 


1 








) 












1 
















1 


1 





1 


1 





1 




1 


















1 






























1 


1 








1 






1 














1 






1 






























1 


1 


1 












1 








\1 












1 





0/ 



(52) 



and so forth. Here all other blanks stand for zeros. We 
omit them just for seeing the sparse matrix structure cor- 
responding to the network geometry of the DSFN. What 
is important here in the above is that the adjacency ma- 
trix for a certain generation of the network can be almost 
3-block diagonal by those for the last generation of the 
network. From this, the fractal nature of the adjacency 
matrices for the DSFN is now very clear. 



VI. EIGENEQUATIONS 

Let us consider the eigenvalue problem for the adja- 
cency matrix in the DSFN. The eigenequation for the 
n-th generation of the network is given by 



A n X n — XX n , 

xi -x n = 1, 



(53) 
(54) 



where A n is the 3™ x 3™ matrix and x n the 3™-dimensional 
vector with its transpose x\ . This reduces to the follow- 
ing eigenvalue problem: 



det[A/„ - An] = 0, 



(55) 



where /„ stands for the 3™ x 3™ unit matrix. Denote by 
D n (X) the determinant in the left hand side such that 

D»(X) = det[XI n - A n ]. (56) 

This can be formally expanded with respect to A as 

D n (X) = A 3 " - aiA 3 "- 1 + a 2 A 3 "- 2 + ■ ■ ■ 



N(X) 



+ (-l) 3 "- 1 a 3 "-iA + 



-l) 3 "a 3 >. 



(57) 



We note the following: Since the A is an eigenvalue of 
the adjacency matrix A n , from the knowledge of linear 
algebra, we are able to derive the following: 



D n (A n ) = A 3 n - ai Al 



a 2 A n 



3" -2 



+ (-l) 3 "-W-iA„ + (-l) 3 "a 3 "/« - 0. (58) 
For example, we easily have the following: 

Di(Ai) = A\ - ai A\ + a 2 A 1 - a 3 h = 0, (59) 



D 2 (A 2 ) = A 9 n2 - aiA\ + a 2 A\ - a^A% + a 4 A 



-a 5 A 2 + a 6 Al - a 7 A 2 + 09/3 = 0, 



(60) 



and so forth. As we will see later, all the terms of even 
number powers of A vanish in the DSFN. This is at- 
tributed to the topology of the network. 



VII. NUMERICAL CALCULATION OF THE 
SPECTRA 

Let us obtain the spectra of the adjacency matrices, 
discussed in the above. We use a computer calculation for 
this purpose. The calculated spectra for the adjacency 
matrix for each generation of the network are shown in 
Fig.2. 

From the numerical results, we find the following very 
important characters of the spectra: (i) The maximum 
eigenvalue, X^\ at the n-th generation of the network 

becomes the second largest eigenvalue, X 2 1+1 ^ , at the (n+ 
l)-th generation of the network. 

(ii) Similarly for the other eigenvalues, all the eigenval- 
ues appeared at the previous generations of the network 
always exist in the eigenvalues appeared at the newly 
generation of the network. 

(iii) The spectrum consists of highly degenerate levels. 
For example, for n = 1 there are three levels (shown by 
red) of A = 0, ±-\/2 with single degeneracy. For n = 2 
there are nine levels (shown by pink) in the spectrum. 
There is only one peak at the center level of A = with 
degeneracy of 3 and the other levels of A = ±(V3 — 
1), ±\/2, ±(\/3 + 1) are all single levels. For n = 3 there 
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FIG. 2: (Color online) The spectrum of adjacency matrix for 
the deterministic scale-free network. The spectrum of adja- 
cency matrix for the network of n-th the generation is shown 
for n — 1 (red); n = 2 (pink); n = 3 (green); n = 4 (blue), 
respectively. 



are 27 levels (shown by green) in the spectrum. There is a 
highest peak at the center level (A = 0) with degeneracy 
of 9. There are two peaks at the levels of A = ±\/2 with 
degeneracy of 3. And the other 12 levels are all single 
levels. For n = 4 there are 81 levels (shown by blue) in 
the spectrum. There is a highest peak at the center level 
(A = 0) with degeneracy of 27. There are two peaks at 
the levels of A = ±\/2 with degeneracy of 9. There are 
four peaks at the levels of A = ±\/2 with degeneracy of 
3. And the other 24 levels are all single levels, and so 
forth. 

In order to study this nature further, we have done 
numerical calculations for the spectra of the adjacency 
matrices with the sizes of 3 5 = 243, 3 6 = 729, 3 7 = 
2187, 3 8 = 6561, respectively. From this we have con- 
firmed ourselves that this nature is numerically exact at 
any generation of the network up to n = 8. We show this 
in Fig.3 for the cases of 3 5 = 243 (green), 6 6 = 729(blue), 
3 7 = 2187(red). 

The above nature is remarkable. It enables us to cal- 
culate the exact sequence of the degeneracies in the spec- 
trum at any generation of the network, apart from find- 
ing the exact eigenvalues. By counting the numbers of 
the levels and its degeneracies, we find the following very 
important result. 

(iv) Let us denote fij the degeneracy of the j-th peak 
in the spectrum for the network of the n-th generation. 
Denote by Qj the number of the eigenvalues having the 
same degeneracy of fij. We find 



fij = 3 J for j = 0, 1, 



_ . 2»- 1 -J for j = 1, 2, • - • ,n-l 
^ J 13- 2"- 1 for j = 



(61) 



(62) 



From this we can check whether or not the above for- 
mula is correct. For this purpose, we just reproduce the 
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2.5 ■ 
2.0 ■ 
1.5 ■ 
1.0 ■ 

0.5 ■ 



-10 




10 




-10 _. . , 10 

Eigenvalue a 



FIG. 3: (Color online) The spectrum of adjacency matrix for 
the deterministic scale-free network. The spectrum of adja- 
cency matrix for the network of n-th the generation is shown 
for n = 5 (green); n — 6 (blue); n — 7 (red), respectively. 



total number of eigenvalues of the adjacency matrix as 
follows. 

n— l 

HjQj = ^oQo + MlQl + ' ' ' + Mn— lQn— 1 

J=0 



1 ■ 2 n_1 • 3 + 3 n_1 • 1 + 3 n_2 ■ 2 • 



= 2 n - 1 -3 + 3(3 n " 2 -l + 3"" 3 -2 



3-2 



n-2 



-3-2 r ' 



3 + 3 



m— 1 on— 1 



3-2 



3". 



1 • 2"~ 2 ) 



(63) 



This is nothing but the total number of eigenvalues in 
the spectrum. Hence, the formula is proved. 

Here we would like to point out the eigenvectors (i.e., 
the states) corresponding to the eigenvalues as well. 
Then we find the following very interesting character: 

(v) The eigenvectors (i.e., the states) are very local- 
ized. The states of A = first appeared in the network of 



the first generation are localized only on the pairs of leaf 
nodes such as the pair of 2 and 3, the pair of 5 and 6, and 
the pair of 8 and 9, etc. The states of A = ±\/2 first ap- 
peared in the network of the first generation are localized 
only on the smallest triples centered at the root hub such 
as the triple of 1, 2, and 3, the triple of 4, 5, and 6, and the 
triple of 7, 8, and 9, etc. The states of A = ±(vo± 1) 
first appeared in the network of the second generation 
are localized within the subnetworks having the size of 
that of the second generation, such as the network of 1-9, 



etc. The states of A = ±y 2(3 ± y/3 ± y/ll ± 6%/3) first 
appeared in the network of the third generation are lo- 
calized within the subnetworks having the size of that of 
the second generation, such as the network of 1-27, etc. 
and so forth. And the number of such subnetworks gives 
its degeneracy. 

Thus, we find that the larger the eigenvalue the larger 
the extent of the eigenvector. Hence, we would like to 
conclude that the eigenvector with the maximum eigen- 
value (say, Ai) is most delocalized while the eigenvalues 
of A = are extremely localized at the leaf nodes in 
the network. This nature explains the meaning of the 
spectrum shown in Fig.s 2 and 3. Therefore, since we 
see similar spectra in most SFN models, we may expect 
that the same holds true for every SFN. This is a very 
interesting point in the study of SFNs. 



VIII. SIMPLE OBSERVATIONS FOR D n {\) 

We now show why so. We are going to treat analyti- 
cally the above result in the previous section. To do this, 
let us consider some interesting nature of D n (X). Let us 
first consider Di(X). 



Di(A) 



A -l -: 
-1 A 
-1 A 



(64) 



We can calculate this and convert it as follows: 





A 


-1 


-1 




A -1 -1 


Di(A) = 




1 A 







-1 A 




1 


A 




-A A 


A -2 - 


-1 




A 

-1 






-1 A 





= A 


-2 
A 


= A(A 2 -2). 





A 







Obviously, this provides three eigenvalues of 

A = 0,±V2. 
Therefore, let us denote it as 

^(A) = A(A 2 - 2) = A/i(A), 

where 

A(A) = A 2 -2. 



(65) 

(66) 
(67) 
(68) 
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Let us next consider D 2 (X). It is given by 

A -1 -1 -1 -1 -1 -1 
-1 A 
-1 A 

A -1 -1 

D 2 (X) = -1 -1 A 

-1 -1 A 

A -1 -1 

-1 -1 A 

-1 -1 A 

(69) 

where all blanks stand for zeros. Doing the same proce- 
dure, we find 



D 2 {\) = A 3 



A -2 -2 -2 

-1 A 

A -2 

A -2 

-1 -1 A 

-1 -1 A 



= A 3 (A 2 -2)(A 4 -8A 2 +4) 

= A 3 (A 2 - 2)(A 2 - 2A - 2)(A 2 + 2A - 2). (70) 
Hence, we obtain 

£ 2 (A) = A 3 /i(A)/ 2 (A), (71) 

where 

/ 2 (A) = A 4 - 8A 2 + 4 = (A 2 - 2A - 2)(A 2 + 2A - 2) 



5 2 (A)/i 2 (A). 



(72) 



In the same way, we find the following for the determi- 
nants -D 3 (A) and D^X): 

D 3 (X) = A 9 (A 2 - 2) 3 (A 4 - 8A 2 + 4) 



(A 8 - 24A 6 + 104A 4 - 96A 2 + 16) 



A 9 {/i(A)} 3 / 2 (A)/ 3 (A), 



D 4 (A) = A 27 (A 2 - 2) 9 (A 4 - 8A 2 + 4) 3 



x (A 8 - 24A 6 + 104A 4 - 96A 2 + 16) 



x(A 16 - 64A 14 + 1104A 12 - 742A 10 + 22112A 8 
-29696A 6 + 17664A 4 - 4096A 2 + 256) 



(73) 



and so forth. Thus, we obtain a series of irreducible poly- 
nomials as follows: 

/!(A) = A 2 -2, 



/ 2 (A) = A 4 - 8A 2 + 4 = (A 2 - 2A - 2)(A 2 + 2A - 2), 



/ 3 (A) = A 8 - 24A 6 + 104A 4 - 96A 2 + 16, 



/ 4 (A) = A 16 - 64A 14 + 1104A 12 - 742A iU + 22112A 



12 v ^O\10 



-29696A 6 + 17664A 4 - 4096A 2 + 256 



(A 8 -4A 7 -24A 6 +40A 5 + 104A 4 -80A 3 -96A 2 + 32A+16) 



x (A 8 +4A 7 - 24A 6 -40A 5 + 104A 4 +80A 3 - 96A 2 -32A+ 16) . 

(75) 

Using these polynomials, we are able to represent the 
sequence of the determinants as follows: 

£>i(A) = A/i(A), 



D 2 (X) = A 3 A(A)/ 2 (A), 
A3(A) = A 9 {/ 1 (A)} 3 / 2 (A)/ 3 (A), 

D 4 (X) = A 27 {/ 1 (A)} 9 {/ 2 (A)} 3 / 3 (A).A(A), (76) 

etc., where / n (A) is an even function of 2™-th order poly- 
nomial of A such that / n (A) = /„(— A). This symmetry 
in the eigenvalues is due to the bipartite nature of the 
DSFN. 

Therefore, we expect that we can carry out such a fac- 
torization of the determinant at any generation of the 
network. Hence, we are led to state the following conjec- 
tures. 

Conjecture 1 For the polynomial with even suffix of 
n = 2m, it is always factorized as 

/ 2m (A) =.g m (A)/i m (A), (77) 
where g m (X) and h m (\) are m-th order polynomials of X. 

Conjecture 2 

A l (A) = A 3 "- 1 {/ 1 (A)} 3 ^ 2 {/ 2 (A)} 3 "" 3 



A 27 {/i(A)} 9 {/ 2 (A)} 3 / 3 (A)/ 4 (A), (74) 



{/„- 2 (A)} 3 /„_ 1 (A)/ n (A). 



(78) 
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The meaning of the second conjecture is now very clear. 
(1) The spectrum is symmetrical around the center level 
of A = (see Fig. 2 and Fig. 3). This means that if A 
is an eigenvalue of the adjacency matrix, then so if for 
—A. (2) The powers in Eq.(78) represent the numbers 
of the degeneracies of the eigenvalues. For example, the 
A = level consists of 3™ _1 (i.e., the degeneracy of 3™" 1 ) 
and this corresponds to the highest peak at the center of 
A = in the spectrum. The A = ±V2 levels consist of 
3™~ 2 (i.e., the degeneracy of 3™~ 2 ) and these correspond 
to the second highest peaks at the center of A = ±V2 in 
the spectrum, and so forth. This proves the formulas of 
Eq.s (61) and (62), previously found. 

The validity of the conjectures is also supported by 
our numerical calculations as mentioned before. But the 
exact proofs have not been made yet, however. 

IX. THE ROOTS OF THE IRREDUCIBLE 
POLYNOMIALS 

We now consider zeros (i.e., roots) of the irreducible 
polynomials of / n (A) for k = 1, 2, • • • ,n. As studied be- 
fore, the order of / n (A) is 2" and it is a function of A 2 . 
Let us denote as y = A 2 . Then, /„(A) becomes a function 
of y such as f n (y)- Now we have the following: 

h{v) = y-% 

My) = y 1 - 8y + 4, 
f 3 (y) = y 4 - 24y 3 + 104y 3 - 96y + 16, 
h(y) =V 8 - 64/ + 1104y 6 - 742y 5 + 22112y 4 



of order up ton- 1, i.e., /i(A), / 2 (A), . . . , /„_i(A). There- 
fore, there are 1 + 2 + • • • + 2 n_1 = 2" — 1 eigenvalues 
already in the spectrum. Thus, the number of newly 
appearing eigenvalues is exactly one more than that of 
the previously existed eigenvalues. Hence, by the knowl- 
edge of algebra, the eigenvalues of the network of the n- 
th generation sandwiches the previously existed eigenval- 
ues such that the maximum eigenvalue may exceed that 
of the previous generation. As we iterate the network, 
n becomes large indefinitely. Therefore, the maximum 
eigenvalue develops further. 

To understand this point further, we have numer- 
ically investigated the ratio R(n) between the maxi- 
mum eigenvalue of A maa; (n) at the n-th generation and 
that of X ma x(n — 1) at the (n — l)-th generation up to 
n = 8. Denote the ratio by R{n) = \ m ax(n) / A max (n— 1). 
The result is as follows: i?(l) = 1.932 R(2) = 
1.583- •• ,R(3) = 1.486- •• ,i?(4) = 1.447- •• ,R(5) = 
1.430 • • • , R(6) = 1.422 • • • , R(7) = 1.418 • • • . From this 
we see that as n becomes large the ratio tends to the 
number \/2 = 1.4142 • • • . Thus, we are led to the follow- 
ing conjecture: 

Conjecture 3 As n — > oo, 

X max -> 2™/ 2 . (80) 

This conjecture can be proved from the nature of the 
series of the irreducible polynomials in Eq.(75). Let us 
first consider the case of n = even. As we have discussed 
Conjecture 1, if n = even, then the polynomial can be 
factorized as 

V -ai{n)y + a 2 (n)y H 

= (A 2 ^ 1 -/3 1 (n)A 2( "- 1, - 1 + -..) 



-29696y 3 + 17664y 2 - 4096y + 256, (79) 
and so forth. 

Since f n (y) is now a 2 7l -th order polynomial of y, it has 
to consist of 2™ zeros. Then the meaning of irreducibility 
is the following: f n (y) (k — 1,2, ••• ,n) does not share 
common roots with the other generations of the poly- 
nomials. This can be proved by the Sturm theorem for 
polynomials . |32j . 

From the knowledge of algebra, if there is a series of the 
irreducible polynomials, then the roots of the polynomial 
of order n — 1 always exist in the intervals between the 
roots of the polynomial of order n. Therefore, the max- 
imum root of the polynomial of order n exceeds that of 
order n — 1. In our problem, the irreducible polynomial 
/„(A) is of order 2™ and it gives the newly appearing 
eigenvalues. Therefore, there are 2™ eigenvalues in the 
network of the n-th generation. 

On the other hand, the previously appeared eigenval- 
ues are given as the roots of the irreducible polynomials 



x (X 2 ^ 1 +p 1 (n)\ 2ln - 1) - 1 + ■■■)= 0. (81) 

where 

I3 l {n) = 2 n / 2 . (82) 
Therefore, the maximum eigenvalue is given by 

A 2 "" 1 - /3 1 (n)A 2< "" ) - 1 ---- = 0. (83) 
Dividing the above polynomial by A 2< >_1 yields 

A = /5i(n)+oQ+-. (84) 

Since \ m ax ~ * co as n — > oo, we can obtain the approx- 
imation of the maximum eigenvalue by a perturbation 
method. Hence, for n — even, we obtain 

A ma ,H = 2"/ 2 -|-0(2( n - 1 )/ 2 ) + -- - , (85) 
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such that X max (n) > 2™/ 2 . Similarly, if n = odd, then 
we can apply the same perturbational argument to the 
irreducible polynomial, then we have 



y = a 1 (n)-0[ - ) +••■ 



(86) 



where a\(n) = 3-2™ and y = A 2 . Hence, we obtain 

X max (n) = V3- 2"/ 2 - O + • • • , (87) 

such that X ma x(n) < V3 • 2™ / ' 2 . This argument sup- 
ports the conjecture. Thus, we expect that as n — > oo 
A ma » - 2"/ 2 . 

In this way, the solution of the series of the irreducible 
polynomials is very crucial in finding the spectrum of the 
DSFN. This point is supported by our numerical calcu- 
lations before. 

X. HIDDEN SYMMETRY IN THE MODEL 

Let us consider a particular nature of the DSFN. There 
is a hidden symmetry in the adjacency matrix, A n . To 
see this point let us consider the case of n = 2. As was 
discussed before, the adjacency matrix for this case is 
given by Eq.(52) and the eigenvalue problem is 
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1 
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1\ 


( X1 \ 
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xs 




x 3 
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XA 




X4 


1 
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X5 


= X 


X5 


1 






1 














x 6 




x 6 
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x 7 




x 7 


1 












1 








x s 




x s 














1 


0/ 


\x 9 J 




\Xg J 



■ (88) 



This expression depends on the choice of the ar- 
rangement of components of the eigenvector X 2 = 
{x\, X2, X3, Xi, X5, xq, x 7 , xs, xqY, where 4 means its trans- 
pose. Reminding ourselves of the bipartite nature of the 
DSFN, we can rearrange them as follows: 



(xi 1 x 4l x 7 ,x 2 ,x 5l x 8l x 3l x e ,x 9 ) t 



(89) 



This means that we align the vector-components num- 
bering in mod3. Then we find the following: 









1 
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1\ 


f xl \ 




( Xl \ 
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X A 




X4 
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x 7 




x 7 
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X2 




X2 
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1 

















X5 


= X 


X5 
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x s 




xs 
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xs 




X3 
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x e 




x e 


\1 





1 












/ 


V x 9 J 




\ x 9 J 



(90) 



where all blanks stand for zeros. Therefore, we can write 
it as 



3 



x3 



M. 



M 2 6x6 



A| I 

v 



(91) 



where u and v are the 3 and 6-dimensional vectors, and 
3X 3 and 6x6 are zero matrices, and 



Ml 



M2 = 



111111 
10 10 
1 1 



(I 0\ 

1 1 

1 1 

1 

1 1 

\1 1/ 



(92) 



(93) 



Here, A 2 = PA 2 P~ 1 and Y 2 = {u,vf = PX 2 , respec- 
tively. This is just a interchange of the vector compo- 
nents. Therefore, if there is no confusion between A 2 
and A 2l we can simply write A 2 as A2. So, we identify 
them as the original adjacency matrix A2, below. 
From Eq.(91), we leads to the following: 



M\v = Xu, M 2 u = Xv. 
From this we find by simple algebra 

X 2 u, M 2 m1v = X 2 v, 



MlM 2 u 



(94) 



(95) 



where M\M2 is the 3-dimensional matrix and M2M\ the 
6-dimensional matrix. These are given by 



M\M2 



111111 
10 10 
1 1 



(\ 0\ 

1 1 

1 1 

1 

1 1 

\1 1/ 



m\M2 = 









(i 






/I 


0\ 


1 


1 


1 


1 


1 





1 


1 


Vi 


I J 




(96) 



111111 
10 10 
1 1 



/I 


1 


1 


1 


1 




1 


2 


1 


1 


2 


1 


1 


1 


2 


1 


1 


2 


1 


1 


1 


1 


1 


1 


1 


2 


1 


1 


2 


1 


\1 


1 


2 


1 


1 


2/ 



(97) 
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Let us solve the eigenvalue problem of Eq.(96). In this 
case, the first eigenvalue problem, M\M 2 u = A 2 u, gives 
the determinant 



My) = 



y -6 -2 -2 
-2 y-2 
-2 y-2 



-0, 



(98) 



where y = A 2 . By simple algebra, we show the following: 



My) 



y-6 -2 -2 
-2 y-2 
-2 y-2 



y-4 -2 -2 
y-2 
-2 y-2 



y-4 -2 
-2 y-2 



(y-2) 

= (y-2) {(y - 2)(y - 4) - 4} = (y-2)(y 2 -8y+4). (99) 

Hence, the determinant Ds{y) gives the three eigenvalues 
of y = 2, \/3 ± 1, as expected. 

Similarly, the second eigenvalue problem, A^M^tT = 
X 2 v, gives the determinant 



My) = 



y-l -1 -1 -1 -1 -1 

-1 y-2 -1 -1 -2 -1 

-1 -1 y-2 -1 -1 -2 

-1 -1 -1 y-l -1 -1 

-1 -2 -1 -1 y-2 -1 

-1 -1 -2 -1 -1 y-2 



= 0. 



(100) 



Subtracting the first column by the fourth column, the 
second column by the fifth column, and the third column 
by the sixth column, respectively, we obtain 



My) 





y-l 


-i -i -i 


1 


-1 




-l 


y-2 -1 -1 


2 


-1 




-l 


-1 y-2 -1 


1 


-2 




-l 


-1 -1 y-l - 


1 


-1 




-l 


-2 -1 -1 y- 


2 


-1 




-l 


-1 -2 -1 


1 


y -2 


y 





-1 -1 -1 









y 


-1 -2 -1 












y -1 -1 -2 






-y 


y-l -1 -1 




( 





-y 


-1 y-2 -1 












-y-l -1 y-2 







(101) 



And as the next step add the first row to the fourth row, 
the second row to the fifth row, and the third row to the 
sixth row, we obtain 



My) = 



-1 
-1 
-1 



y 
y 
y 
y-2 
-2 
-2 



-1 
-2 
-1 
-2 
y-4 
-2 



-1 
-1 
-2 
-2 
-2 
y-4 



= y 



y-2 -2 -2 
-2 y-4 -2 
-2 -2 y-4 



y-2 -2 
-4 y-6 
-2 2-y y-4 



y-2 -2 
-2 y-2 -2 
-2 2-y y-4 



y 3 (y-2) 



y-2 -2 
-4 y-6 



y 3 (y-2)(y 2 -8y + 4). 



(102) 



Hence, the determinant D e (y) gives the six eigenvalues 
of y = (triple), 2, \/3 ± 1, as expected. We would 
like to note that the factors seen in Eq.s (99) and (102) 
are nothing but the irreducible functions /i(y) and /2(y), 
found before. 

Let us consider some interesting character in Eq.(91). 
Since A 2 X 2 = XX 2 , we find that 



(A 2 ) 2 X 2 = X 2 X 2 , 



(103) 



where 



3x3 Ml 
M 2 6x6 



3x3 Ml 
M 2 6x6 







M\M 2 





m 2 m1 



(104) 



Hence, the square (A 2 ) 2 of the original adjacency ma- 
trix A 2 for the DSFN can be block-diagonal. Since from 
Eq.(103) we have {{A 2 ) 2 - X 2 }X 2 = 0, this eigenvalue 
problem provides the following: 

det [{A 2 f - \ 2 I 2 ] = dot [A 2 - XI 2 ] dot [A 2 + \I 2 ] 



= I det [A 2 - XI 2 ] 



det 



m\m 2 - X 2 h 



det 



M 2 M\ - X 2 h 



(105) 



Therefore, the original determinant D 2 (X) of Eq.(69) is 
given by 



D 2 (X)= det M\M 2 - X 2 I 2 det M 2 M\ — X 2 I 2 



1/2 



1/2 



= (My)Mv)) 1/2 - (ioe) 

Substituting Eq.s (99) and (102) into Eq.(106), we obtain 
D 2 (X) = y 3 / 2 (y - 2)(y 2 - 8y + 4) = A 3 / 1 (A)/ 2 (A). (107) 
This reproduces Eq.(76). 
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Now, we are able to generalize the above procedure to 
the adjacency matrix A n for the DSFN at any generation. 
In this case the Eq. (91) is generalized to 



Ouu M\ 
M n VV 



= A 



This then yields 



M^v — Aw, M n u = Xv. 



(108) 



(109) 



Here u and v are the u- and v-dimensional vectors, MU 
and VV the u x u and v x v zero matrices, and M\ and 
M n the u x v and v x u matrices, respectively, where 
u = 3 n_1 and v — 3" — 3 n_1 . Then, in the same way we 
have the following: 



D n (X) = |det [MlM n - A 2 /„] | V2 |det [M n M+ - A 2 / n ] 



1 1/2 



D 3 „-i(y)-D 3 „_3„-i(y) 



1/2 



(110) 



where D 3 n-i(y) and L> 3 „_ 3 „-i (y) [= y 3 D 3 „-i(y)] are 
factorized in terms of the irreducible polynomials of 
Eq.(79). Therefore, Eq.(107) may reduce to the form 
of Eq.(75). Thus, this approach may provide a hint to 
prove the conjecture 2. However, the proof is out of scope 
of this paper since it is numerically supported as in the 
section VII. 



We would like to mention that similar ap pro aches 

r pho t 

topological localization problem 35], Fermion number 



j pro acnes 

have been applied to the amor pho us systems]^ 0], 

l r 

fractionalization[36j, and the Hubbard model 37] . 



XI. ZERO MODES AND INDEX THEOREM IN 
THE DETERMINISTIC SCALE-FREE 
NETWORK 

Let us consider zero modes (i.e., the eigenvectors hav- 
ing A — 0) in the spectrum. The zero modes are given 
by substituting A = into Eq.(107) such that 



M}.v = a, M n u = 0. 



(Ill) 



The number of zero modes given by M n u — (or equiv- 
alently, the number of zero modes given by the matrix 
M^Mn) is the dimension of the null space of M n . This 
is sometimes called the dimension of the kernel of M n , 
that is simply written as dimjker M n ], And the number 
of zero modes given by M^v — (or equivalently, the 
number of zero modes given by the matrix M n M^) is 
the dimension of the null space of Aft. This is some- 
times called the dimension of the kernel of Aft, simply 
written as dim[kerAft] = dim [co ker M„ ] |3q . Therefore, 
in the above example of n = 2, dim [ker M 2] — 3 and 
dimpcerMj] = 0. This would give a relation: 

Ind^Mz) = the number of zero modes 



dim[kerM 2 ] - dim[ker Af 2 T ] = 3 



(112) 



This quantity Ind(M n ) is called the index of M n . And 
the relation that the number of zero modes coincides with 
the difference between dim[kerM n ] and dim[kerAf^] is 
called the index theorem\jtffy . 

This can be generalized to the adjacency matrix A n 
for the DSFN in the n-th generation such that we have 

Ind(M n ) = dim[kerM„] - dim[ker M*]. (113) 

As we have obtained Eq.(78), for this case we obtain 

dim[kerM„] = 3"~\ dim [ker A/,^] = 0. (114) 

Hence, we derive the following index theorem for the 
DSFN: 

Ind(M n ) = dim[ker M n ] - dim[ker Aft] = 3"" 1 . (115) 

As we have discussed before, this number is just the num- 
ber of A = 0, localized at the smallest leaf nodes. 

The proof of Eq.(114) is quite simple. For the DSFN, 
the dimension of the matrix M^Mn is 3™ _1 , which is 
the number of the hub nodes. And the dimension of the 
matrix M n M t is 3™ — 3" _1 , which is the number of the 
rim nodes. Then, we always find that there are no zero 
modes in AftAf„ in our DSFN. That is, dim[kerMt] = 
0. Therefore, the null space of M„Aft is given by the 
difference between the dimension of the matrix M n M t 
and that of the matrix AftM„. This is 3"-3 Tl " 1 -3™^ 1 = 
3" _1 . Hence, dim[ker MJ = 3™" 1 . 



XII. THE NATURE OF THE MAXIMUM 
EIGENVALUE 

Let us consider the nature of the maximum eigenvalue. 
As we studied before, the matrix M^M n absolutely de- 
termines the spectrum of A > 0, while the matrix Af n Aft 
determines the zero modes of A = as well as the spec- 
trum. Therefore, in order to determine the spectrum we 
need study the matrix M\M n . 

Let us see some examples for this matrix. The matrix 
for n = 2 is given in Eq.(96). We also give the matrix for 
n = 3 as 



M 3 T M 3 



(U 

2 
2 
4 
2 
2 
4 
2 

V 2 



2 2 4 2 2 4 2 2\ 
20000000 

2 2 



2 

6 

2 

2 




2 
2 



6 
2 
2 



2 2 
2 
2/ 



(116) 



Here, we find that the diagonal elements of the matrix 
AfjA/ 3 are just the maximum degrees of nodes (i.e., the 
numbers of the most connected links) in the networks up 
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to the third generation. The maximum diagonal element 
is the maximum degree of nodes, which is the degree of 
the root. Hence, it is 14 for n = 3 [Obviously, it is 6 for 
n = 2 as in Eq. (96)]. 

We find that this is always true for the network of 
the n-th generation. We find that the diagonal elements 
of the matrix M\M n are just the maximum degrees of 
nodes for the networks up to the n-th generation. Hence 
the maximum diagonal element is the maximum degrees 
of node, which is the degree of the root. It is k max = 
2(2™ - 1) from Eq.(7). 

Let us first derive the lower bound of the maximum 
eigenvalue, Amax- In mathematics, we have the following 
theorem |39|: 

Theorem 1 Denote by H a non-negative definite ma- 
trix. Denote by v a positive definite symmetric ma- 
trix. Define a matrix A by A = H + V . Suppose that 
the eigenequations H\ipi) = hi\ipi) and A\ipi) = a^ipi) 
(i = 1,2, ■ ■ ■ ,1) are known such that the eigenvalues sat- 
isfy hi > h% > • • • > hi and a± > 02 > ■ • • > a; . And 
assume that (tpi\V\ipi) > 0. Then, the following holds 
true 
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FIG. 4: (Color online) The growth of the maximum eigen- 
value, Xmax and its lower and upper bounds in the DSFN. 
Here Xmax (circle), the lower (triangle), and upper (square) 
bounds are shown for the DSFNs of the n-th generations up 
to n = 8, respectively. 



h s < a s for s = 1, 2, 



■ I- 



(117) 



We omit the proof here, since it is given in the 
literature [39^. 

Let us apply this theorem to the matrix M\M n . We 
now denote the matrix M\M„ by A in the theorem. De- 
note by H the matrix having the diagonal matrix compo- 
nents of M\M n and denote by v the matrix having the 
off-diagonal matrix components of M^M n . Therefore, 
M\M n satisfies the conditions for A in the theorem. The 
eigenvalue of the matrix A is now a* = Af and I — 3 n_1 . 
Therefore, we can prove the following theorem for the 
DSFN. 



Theorem 2 



h s <X 2 s for s = 1,2, 



Since the maximum eigenvalue of H is now k max , we have 
hi = k max . Thus, we finally end up with the following 
theorem for the DSFN. 



Theorem 3 



\fk — \/2 n + 1 — 2 < X 



(119) 



We have also investigated this theorem numerically in 
Fig. 4. It shows that the theorem is valid for the DSFNs 
of the generations up to n = 8. Thus, this supports the 
theorem. 

Let us next derive the upper bound of the maximum 
eigenvalue, X max . To do so, let us first consider a par- 
ticular property of the matrix, M n . As is seen from the 
matrix such as Eq.(93), the i-th column vector of M n are 
a vector whose components are or 1 such that the total 



number of 1 in the column counts the order ki of the i-th 
hub node. Denote this vector by ki. Then we find 



kj, ' ki ki 



(120) 



where by definition k\ > k2 > • • • > k%n , since k\ is the 
order of the root (i.e., the most connected hub). Then, 
the matrix can be represented as 



M n = (ki,k 4 ,--- ,fc 3n _ 2 ), 



(121) 



where the suffices run over all the hub (i.e., root) nodes. 
Now, we can represent M\M n in general in the following: 



(US) MlM n = 



( k\ ■ k\ k\ ■ ki 

k\ ■ ki k^ ■ ki 



fel-fe3"-2 \ 

ki ■ fe 3 n_ 2 



\ kl ■ k 3 n_ 2 k 2 -k 3 n_ 2 ■■■ fc 3 "-2 • fc3"-2 / 

(122) 

which is a 3™ 1 x 3" 1 symmetric matrix. Therefore, the 
trace of this matrix is the total number of links such that 

tr (M^Mn) = h-ki= k i 



i£root 



i^root 



J2 kP{k) = L{n). 



(123) 



kGroot 



We now need the following theorem known as the 
Perron- Frobenius theorem in linear algebra |3lll39| . 

Theorem 4 Suppose that an n x n symmetric matrix A 
has all non-negative entries aij > 0. Then this satisfies 
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an eigenequation A\tpi) = ai\tpi). For any positive con- 
stants Ci, C2, . . . , c n , the maximum eigenvalue a max (A) 
satisfies 



(124) 



a max (A) < max < — 

Ki<n I ^— ' C 



3=1 



We omit the proof here since this is very well-known |3 ll 

EJ 

Let us apply this theorem to our problem. We now 
regard A as M\M n such that = fcj • kj and a, = Af . 
Then, from the theorem we find the following: 



AL(MlM K ) < max { ^ > ■ ' ' - >r > ' 

lEroot | * — ' Ci 



From this we also find 



X 2 max {MlM n ) < max <j £ 

j (Eroot 



< max { y 

ieroot I ^— ' 



< max 

Ki<» 



j£root 



(126) 



Since |cj/cj| < 1, then we have 



*Lx(^/„) < max I 

i£root | f rf 
j£root 



(127) 



Hence, we have the following theorem for the DSFN. 
Theorem 5 

\ 2 max (MlM n ) < h-kj = b. (128) 

j€root 

Since £ ieroot & ' ^ ^ Enroot ^7 ' % = L ( n ) together 
with Eq.(6), we end up with the following theorem for 
the DSFN. 

Theorem 6 



X max (MlM n ) < y/Lfij = v/2(3« - 2"). (129) 

From theorem 3 and theorem 6, we finally derive the 
following theorem for the DSFN. 

Theorem 7 The maximum eigenvalue of the adjacency 
matrix of the DSFN is bounded as 



< Amai < b < 



(130) 



where a = y / 2{2 n - I), c = ^2(3™ - 2 n ), and b = 

J2 3 (=root ^1 ' kj. 



Although there is no name on the quantity b, this is a 
much better bound than c, the square root of the total 
number of links. We note that the above theorem is 
consistent with Conjecture 3 that asserts X max — * 2™/ 2 
asn^ oo. 

We finally make a comment on the Chung, Lu and 
Vu's theorem |3l| - As introduced in Introduction, since 
the DSFN has the exponent of 7 = In 3/ In 2 < 2.5, this 
system may belong to the (C2) case. Therefore, if we 
apply their theorem to our system of the DSFN, then it 
leads us to the following: 



ykmax — X max — d : 



(131) 



where d oc (4/3)" = (1.33- ••)" horn Eq.(49). On the 
other hand, our theorem provides us the lower bound 
a oc 2™/ 2 = (1.414 • • • )™ and the upper bound c cx 3™/ 2 = 
(1.732 ■••)". From this, we find that the Chung, Lu and 
Vu's theorem |3lj| holds true until about n — 12. However, 

beyond n = 12 the role between d and y/k ma x is switched. 
So, their theorem can be violated above n = 12, although 
strictly speaking, their theorem should be applied to the 
case of 2 < 7 < 2.5. We show this behavior using the 
above exact expressions in Fig. 5. As is shown in Fig. 4, 
the numerical value of X max is sandwiched between the 
upper bound of d and the lower bound of \fk max - There- 
fore, it follows the Chung, Lu, and Vu's theorem within 
our calculation up to n = 8. Unfortunately, the numer- 
ical value of X m ax is not available above n = 8 because 
of the ability of our computer power. So, we cannot say 
anything about where it is located above n — 8, so far. 
We may expect that it is still sandwiched between them 
although the role is switched. 

In this way, our rigorous approach provides a concrete 
example for investigating the validity of the Chung, Lu 
and Vu's theorem[3l| . This is an advantage of our theory. 



XIII. CONCLUSIONS 

In conclusion, we have intensively studied the DSFN 
that was first studied by Barabasi, Ravasz, and 
Vicsek.28]. We have first studied the geometry of the 
network and presented the exact numbers of nodes and 
degrees. From these we were able to calculate the ex- 
act scaling exponents for the hub and rim nodes, where 
we have shown that the scaling nature of the hubs is a 
scale-free law but that of the rims is not so but the expo- 
nential. This is obtained only when we have known the 
exact numbers of the nodes and degrees. So, the scaling 
behaviors of the hub and rim nodes are different. In some 
sense we may say that the DSFN is more like a multifrac- 
tal. This nature has not found yet by Barabasi, Ravasz, 
and Vicsekp^. 

Second, we have analytically calculated the exact num- 
ber of the second order average degree d [see Eq.(49)], us- 
ing the numbers of the nodes and degrees. This quantity 
has been frequently used in the criterion for investigating 
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FIG. 5: (Color online) The growth of the square root of the 
total number of links \J L(n), the maximum order of nodes 
y/kmax, and the maximum degree d(n). Clearly, we see that 
\/k m ax is smaller than d(n) but very close to it under n = 12. 
However, y^maT exceeds d(n) about n = 12. 



the bounds of the maximum eigenvalue X max of a SFN as 
Chung, Lu, and Vu|3ll| have been emphasizing. However, 
there has been no example of a rigorous result. So, this 
would be a first example of the analytically calculated 
quantity, d. 

Third, we have numerically calculated the spectra of 
the adjacency matrix A n for the network up to the n = 7- 
th generation. From this we have counted the exact num- 
bers of degeneracies in the spectrum, and have shown 
that the density of states looks like a fractal, where there 
is a large peak at the center of the spectrum of A = [see 
Fig. 2 and Fig. 3]. Such degeneracies are related to the 
degree of the localization of the state (i.e., the eigenvec- 
tor) . As the eigenvalue becomes larger the state becomes 
broader (i.e., delocalized) , where the minimum eigenval- 
ues of A = are localized on the minimum connected rim 
nodes. The envelope of the density of states in the DSFN 
is somehow similar to that of the spectral shape in the 
AB-model0. So, we expect that the origin of localiza- 
tion in the AB-model may be the same or similar to that 
in the DSFN that we have studied here. 

Fourth, we have discussed the nature of the adjacency 
matrix for the network. We have shown that there is a 
recursive structure of the determinant, which induces to 
the complex structure of degeneracy. The sequence of the 
irreducible polynomial seems very important. This would 
be related to a new type of functions. The roots of the 
polynomials and their irreducibility may be a problem of 
the Sturm theorem in algebra. Therefore, this will be an 
interesting problem for mathematicians. 



Fifth, we have shown that there is a hidden symmetry 
in the adjacency matrix. This is the consequence of a 
bipartite structure of the network. From this nature, 
the adjacency matrix is decomposed into the off-diagonal 
type. Then we have proved an index theorem for the 
network for the first time. This theorem has enabled us 
to count the number of zero modes exactly. This type of 
theory is very well-known in the quantum field theory [3tl 
l38[ . It is known as the supersymmetric structure. And 
it has been also known in solid state physics from a long 
time ago [33, E3 HE HE H3- The number of zero modes 
is related to the number of localized states in the system. 
Therefore, we expect that this may be the case in other 
network systems. This will be a very important problem. 

Finally, we have investigated the maximum eigenvalue, 
Amaz in the DSFN. We have shown that X max is bounded 
by the lower and upper bounds that are given analyt- 
ically [see Theorem 7]. These bounds grow as fast as 
exponential as the network grows. Hence, by the the- 
orem, X max grows exponential as well. This might be 
also a first example for the determination of the exact 
bounds. In the standard networks such as the random 
networks |2(, the maximum eigenvalue cannot grow so 
fast as the network grows. And also, as in solid state 
physics, networks in most of physical systems provide 
the so-called energy band that is a spectrum with a fi- 
nite region. This is due to the topology of the finite 
coordination number of atoms in the network of the lat- 
tice structure. So, in order to elucidate the difference 
between the SFNs and other networks the growth of the 
maximum eigenvalue is an important signature. As was 
studied by many authors [J, 011 HE EES}, the behavior 
of Xmax in the AB-model is believed to be proportional 
to y/kmax- Therefore, since k max is propotional to V~N, 
hence, X max oc N 1 / 4 . To see whether or not this is true in 
an arbitrary SFN and to know how general it is, Chung, 
Lu, and Vu[3l| established the theorem. However, as we 
have shown, X max in the DSFN that we have studied in 
this paper seems to escape from the region of their the- 
orem. Therefore, our example will be a counter example 
for their theorem. Whether or not this is true will be an 
interesting problem for further study. 

Thus, we believe that our theory presented in this pa- 
per gives the first rigorous example in the SFN theory, 
where most of all quantities in the network theory are 
analytically obtained. In this context, we would like to 
call the DSFN the exactly solvable SFN model. 
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